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The Principles of Internal Ballistics. 
By J. Pkoudman, University of Liverpool 

(Communicated by H. Lamb, F.E.S. Eeceived June 11, 1921.) 

1. This paper contains the results of an investigation carried out at the 
Eesearch Department, Eoyal Arsenal, Woolwicli. It is published by per- 
mission of the War Office. The principal features are the methods of 
determination of the rates of burning of propellants and the treatment 
of the mathematical equations relating to the motion of the shot in the 
gun. No attempt has been made to improve on the methods of treating the 
distribution of motion and pressure through the gases in the gun, the resist- 
ances to the motion of the shot^ the conduction of heat through the metal 
of the gun, or the recoil of the gun. 

The general theory is sketched in terms of functions which require experi- 
mental evaluation • and then assumed provisional forms are assigned to these 
functions. 

The author wishes to state that he owes the foundation of his knowledge of 
the subject, and also the opportunity of making these investigations, to 
Messrs, F. B. Pidduck and A. D. Crow. 

Gmwml Theoretical Considerations, 

2. We shall assume that the density and pressure of the propellant gases 
determine the temperature, so that we have always an equation of the form 

T = H(A,P), (1) 

where A denotes the density, P the pressure, and T the absolute tem- 
perature of the gases. We shall also assume that the chemical constitution 
of the gases is determined by the density and pressure, the attainment of 
chemical equilibrium being so rapid that the state appropriate to any pair of 
values of A and P is reached before these values have been sensibly departed 
from. We can then treat the gases arising from the burning of a propellant 
as a single body, whose state is defined by A and P, and all absorption or 
emission of heat due to changing chemical constitution may be allowed for in 
the specific heats taken. 

Let E(A, P) 

denote the intrinsic energy per unit mass of the gases when in the state 
given by A and P. By intrinsic energy we mean the maximum amount of 
work which the gases can be made to do by allowing them to expand to zero 
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pressure and to cool to absolute zero of temperature. Let E denote the 
intrinsic energy per unit mass of the gases before they have done any work 
or lost any heat after being generated from the burning propellant. This we 
may call the intrinsic energy per unit mass of the unburnt propellant ; it 
will depend only on the state of the unburnt propellant. The "uncooled 
ipressure-density relation '' will be 

E (A, P) = E. (2) 

If we plot P against A, then the possible states of the gases will be given 
'by all the points of an area, and (2) will be the equation of the upper 
boundary of this area. 

Let dv and a-p denote the specific heats at constant volume and constant 
pressure respectively ; these will both be functions of A and P. If we write 

^=1/A, (3) 

iso that V is the volume per unit mass of the gases, we have, J denoting the 
mechanical equivalent of heat. 
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Van der Waal's theoretical form for the function H (A, P) is 

EH (A, P) = (P + ^A2)(1/A-Z>), (6) 

where E, a, &, are positive constants. If this law applied to the gases in 
question and if a^ and a^, were constants, we should have from (5) 

(Tp — crv = E/J, (7) 

and then the integration of the equations (4) would give 

(7~-l)E(A, P) = (P + 7aA2)(l/A~-&), (8) 

where 7 = o-p/av 

The uncooled density-pressure relation (2) would become 

(P + 7aA2) (1/A-h) = n, (9) 

.where n = (7-~l)E. (10) 
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The equation (9) may be written in the form 

so that if we plot P/A against P, we obtain a curve which has for an 

asymptote 

P/A = &P + n; (12) 

and always lies below this asymptote, except on the axis of P = 0, where the 
two meet. Equation (12) expresses Noble and Abel's pressure-density 
relation. 

When a piece of homogeneous and isotropic propellant is surrounded by 
gas in a uniform state, we assume that combustion takes place by parallel 
layers. If u denotes the rate of hurning down any normal, we shall suppose 
it determined by the state of the gases, so that 

u = v{^, P). (13) 

Let us consider an explosion in which all the constituent pieces of the 
charge are exactly alike, have been ignited all over their surfaces at exactly 
the same instant, and the state of the gases in contact with the unburnt 
charge is continuously uniform. 

Let vq denote the volume and so the total surface of a piece of propellant 
before ignition. At any time, t, during the explosion, let z denote the 
fraction of the charge burnt, and s the total surface of each piece. We 
shall have 

S=: So(l> (Z), (14) 

this defining (f> (z), Charbonnier's form-function, and also 

Vq dz/dt = us, (15) 

which give 

l=f:-^(4 (16) 

In ballistic applications we require the functions 

E-E(A,P), U(A,P), 

for all the values of A and P that occur. The function H(A, P) and the 
specific heats cr^ and <rp may occur as auxiliaries in ballistic investigations, 
but, on the present assumptions, need not enter explicitly into practical 
ballistic calculations. 

Closed Vessel Explosions. 

3. Let us first consider an ideal closed vessel with non-conducting walls, 
and let us suppose that the attainment of mechanical equilibrium is so rapid 
that A and P may be taken as uniform through the vessel at any stage of 
the explosion. 

u 2 
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Let Ai and Pi denote the maximum final values of A and P respectively. The 
relation between A and P will, throughout the explosion, be given by (2) and 
may be found from that between Ai and Pi in different explosions. If S 
denote the density of the unburnt propellant, supposed uniform and constant 
throughout the explosion, we have 

l/A-l/S' ^^^^ 

it being supposed that there is no non-gaseous residue after the explosion. 

We suppose that the record gives P as a function of t ; then (2) gives A 
and (17) ^ as functions of t. Using (16), we obtain to as a function of t, z, A 
or P. It is to be noted, however, that we have only considered here the 
determination of the function U (A, P) for the gases in the uncooled state ; 
that is, only for the upper boundary of the representative area in the A, P 
plane over which the function is required. 

If, after the explosion, we could draw off the heat from the gases so as to 
keep Ai constant and A and P uniform through the vessel, and could take a 
record of the quantity of heat, Q, given off, as a function of P, we should 

JQ = E~-E(Ai, P), (18) 

so that E— E (A, P) could be determined for all desired states. 

Now let us suppose that the release-valve of the vessel is connected with 
that of another vessel whose walls are conducting. Let both vessels be 
charged and the valve kept closed until the charge in the conducting vessel 
has been exploded and the gases allowed to cool. Now let the valve be 
opened, so that some of the cooled gases pass into the non-conducting vessel 
and reach a density Ao and pressure Po round the unburnt charge of 
density Ai. Next let the valve be closed and the charge fired, and suppose 
it burns sufficiently slowly for all the gases in the vessel to be continuously 
in a uniform state. If this state be given by A, P, we shall have 

/ Ao\ l/Ai-l/8 
^-V AJ l/A-1/8' ^^^> 

i-|){E(A, P)-E(Ao, Po)} ={^-^ {E-E(A, P) }. (20) 

When the function E~E (A, P) is known for all the values of A and P in 
question, the equation (20) gives the relation between A and P during the 
second explosion, the gases being now in a partially cooled state. The record 
for this second explosion is supposed to give P as a function of t] by means 
of (20) we obtain A as a function of t, and by means of (19) z as a function 
of t. On again using (16) we can then calculate the rate of burning, u, for 
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states other than uncooled ones. A means is thns indicated for the deter- 
mination of the function U (A, P) for all values of A and P required. 
If we assume 

E(A,P)=^P(H-)' <21) 

which is obtained from (8) by the usual practice of taking a = 0, b = 1/8, 
we obtain from (20) 

P-P n(l/Ao-l/A) .ooN 

» (l/Ao-l/S)(l/A-l/8)' ^ ^ 

as the relation between A and P during the second explosion. 

Cooling, 

4. So far we have imagined ideal conditions. When we come to consider 

explosions in actual vessels, the loss of heat by conduction through the walls 

is a serious factor. It is obvious that with the same density of loading in 

the same vessel the maximum pressure will decrease as the duration of the 

explosion increases, and the effect has been noticed ever since accurate 

experiments were made. It is also clear that for those shapes of pieces of 

propellant which make 

<jb(^)-*0, as 2; -►I, 

there will, at the instant of maximum pressure, be a certain fraction of the 
charge unburnt. For maximum pressure occurs when the rate of gain in 
pressure due to burning is balanced by the rate of loss in pressure due to 
cooling, and for the shapes in question the rate of gain by burning vanislies 
with the fraction of the charge unburnt. Further, for the same rate of loss 
by cooling, the fraction unburnt at maximum pressure will be greater the 
greater the duration of the explosion, or, practically, the greater the size of 
the piece of charge. The fall in maximum pressure from one explosion to 
another with the same density of loading, but with a larger size of piece of 
propellant, is therefore in excess of that representing a loss of energy through 
the walls of the vessel. 

There exists, however, an important cooling phenomenon of a different 
type. This appears to have been first noticed by K. Blochmann* for high 
explosives. He found that for the same density of loading the pressure 
increased as the value of S/C decreased, S denoting the total internal surface 
of the explosion vessel and C its capacity, and he determined the limiting 
value of the pressure as S/O-*0. 

For propellants, the diminution of maximum pressure due to cooling had 

■^ * Dingler's Polytechnisches Journal, vol. 318, p. 234 (1903) ; the author owes this 
and many other references to Mr. F. B. Pidduck. 
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been elucidated in the Department. As was to be expected, it had been 
found that, for the same vessel and value of Ai, the value of Pi increased as 
the duration of the explosion decreased. For each value of S/C and Ai 
taken, the limiting value of Pi, as the duration tended to zero, was estimated. 
It was then found that, for the same value of Ai, these limiting values of Pi 
increased as the values of S/C decreased, this being of the same nature as 
Blochmann's result. For each value of Ai taken, second limiting values of 
Pi as S/C-> were found. These second limiting values of Pi were taken to 
represent the results of an ideal explosion in which no cooling takes place, 
and from them the uncooled pressure-density relation (2) was constructed. 
Further experiments of the kind are urgently required, great attention being 
paid to the temperature conditions of the vessel and propellant, and to the 
regularity of size and shape in the pieces of propellant. 

Apparent Bates of Burning. 

5. During an actual explosion, we do not know the true relation between 
A and P, and so z cannot be found accurately as a function of t for the 
purpose of calculating the rate of burning by (16). 

The method indicated at the beginning of § 3, when applied as if no 

cooling took place, is practically that established by Vieille,* and hitherto, in 

no available publication on rates of burning, have cooling effects been 

eliminated or allowed for. Previous to Vieille's work, writers on ballistics 

had usually taken 

iv = %'V% (23) 

u' and a being constants. For a different writers had used 

1 -l 3. 2. 1 

^j 4; 2? 5^ 3> -^^ 

Vieille decided on the same law, as an approximation, for the propellants 
with w^hich he experimented, and took for a 

2 5 1 
3"' 9 ^ 5 ? 

for nitrocellulose, cordite, and ballistite respectively. 

For all the French propellants, Charbonnier, Gossot and Liouville claim 
the validity of (23), but whereas Charbonnierf maintains that a = 1, G-ossot 
and Liouville:j: insist that « = f . For cordite M.D., heated to 80° F. before 
ignition, Mansell§ found the relation 

It z=z ii^-\-itiV, (24) 

uq and ui being constants. 

* * M6morial des Poudres et SalpStres,' vol. 6 (1893). 

t See, e,g., ' Balistique Interieure,' p. 76 (Paris, 1908). 

X See ' M<3morial des Poudres et Salpetres,' vol. 13 (1906). 

§ * Phil. Trans,,' A, vol. 207, pp. 243-262 (1907). 
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Both (23) and (24) make to independent of A, and it is impossible to 
obtain -i^ as a function of more than one independent variable by the methods 
that have been used. By these methods, each record leads to a relation 
between u and P, which, when the duration of the explosion is not too 
long, is an approximation to the true relation with A connected to P 
according to (2). 

It is clear, however, that as for the same value of P in different explosions 
we may have different values of A, we may expect the derived relation 
between u and P to vary from the analysis of one record to that of another. 
On the one hand, the relation between A and P cannot always be that 
assumed in the analysis, and, on the other hand, the true rate of burning 
may vary with the value of A for the same value of P. 

A first step in the study of true rates of burning is the elimination of 
cooling effects during the whole of the explosion. This requires an 
elaborate series of experiments, such as those already referred to. When the* 
customary methods are used, we refer to the results as " apparent.'' 

An apparent rate of burning is stated as a function of the pressure only ; 
this function varies with the circumstances of the explosion, such as the size- 
and shape of the vessel, the size of the pieces of propellant, the density of 
the charge, the temperature of the vessel and of the charge. 

Uncooled Bate of Burning. 

6. We shall now show how to apply to the whole of the pressure -time 
records a double limiting process, so as to eliminate cooling effects. 

We still make all the assumptions stated at the beginning of § 3, except 
that relating to cooling. From each record we calculate yvo/so as a fraction 
of 3/ = P/Pi. In the absence of cooling effects, this relation would depend 
only on Ai and be independent of Vq/sq and S/C. On this account, and also 
because it obviates the uncertainty of corresponding times in different 
explosions, it is this relation between yvo/so and y, and not that between 
P and t, which is subjected to the limiting process. 

We suppose that we have pressure-time records for each of m values of 

Vo/so in closed vessels having n values of S/C. For each value of S/0, and 

each of the values 

0-1, 0-2, 0-3, ... 0*9,1, 

of y, yvo/so is plotted against the values of Vq/sq, 

From each set of m points thus obtained, the limiting value of yvo/so as 
Vq/sq-^O is estimated. We are thus enabled to plot the limiting relation 
between yvo/so and y for each value of S/C. 
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Again for each of the values 

0-1, 0-2, 0-3, ... 0-9, 1 

of y, the limiting values of yv^js^ already found are plotted against the 
corresponding values of S/C. From each of the set of n points thus 
obtained, the second limiting value of yvo/so as S/C -> is estimated. 

We are thus able to plot the double limiting relation between yvofso and y, 
and this is taken as corresponding to an explosion in which there are no 
€ooling effects. We now suppose that our symbols refer to this hypo- 
thetical explosion, and use the method sketched at the beginning of § 3. 

From (17) z can be calculated as a function of y, since A is connected 
with P by the " uncooled " relation. From (16) we have 

U:=:V^-^~ (25) 

^s,^{z)dy ^ ^ 

from which u can be calculated as a function of y, and therefore of P. 
This represents the '' uncooled " rate of burning function, i.e., the function (13), 
when A and P are connected by (2). 

Vessel of Variable Capacity.^ 

7. Let unit charge be placed in a vessel with non-conducting walls whose 
capacity is variable, and let z denote the fraction burnt when the capacity 
is V. Then we have 

^-(l~-^)/S' 

still regarding S as uniform and constant. 
The total intrinsic energy in the vessel is 

(1-~^)E + ^E(A, P), 



and we deduce 



~^{E(A, P)-~-E}+P = 0. (26) 

dv 



This is the fundamental differential equation for the circumstances we are 
now considering, and if v be interpreted as the capacity of the vessel per unit 
mass of the charge, it is of general application. It is an equation between v, z 
and P which cannot be integrated without further relation between these 
variables. 

The non-integrability of the equation (26) may be seen from the considera- 
tion that the change in the state of the gases for given changes in capacity 
and amount burnt, must depend on the way these changes have been made, 
and not simply on the extent of the changes. If an expansion takes place 

^" Cf. A. G. Hadcock, *Koy. Soc. Proc,,' A, vol. 94, pp. 484-488 (1917). 
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with a small fraction burnt, a small amount of work is done by the low 
pressure ; if the same expansion takes place with a large fraction burnt, a 
large amount of work is done by the high pressure. In these two instances 
the final intrinsic energies of the gases must be different, though the capacities 
and fractions burnt be the same. 

In the absence of experimental knowledge of the function E (A, P) we are 
^gain forced on the provisional assumptions. Assuming the formula (21) we 
have 

^E(A,P) = -l^p(^-ij, (27) 

l-i = i(,-l=i)-l = l(._l). (28) 

On using (iO) the equation (26) becomes 

or n^ = (^'y--iV^ + 'yP (30) 

Equations of the Gun. 

8. Let C now denote the capacity of the chamber of the gun, A the area of 
bore, and m the travel up to time t ; let m denote the effective mass of the 
shot, c that of the charge, and z the fraction burnt up to time, t; also, let V 
•denote the velocity of the shot and P the pressure of the gases at the time t. 
Then we have 

v = (G-hAx)/c, (31) 

Y = dx/dt, (32) 

and m dY/dt = (1 -0) A (P -K), (33) 

where E and are inserted to allow for resistance and the angular accelera- 
tion of the shot. The equation (26) yields 

J- ~ z {E (A, P) -M} + P = 0, (34) 

where A = r^ , (35) 

(G + A^)/c-(l-.^)/8' ^ ^ 

and h is inserted to allow for loss of heat through the walls of the gun. 

The equations (32), (33), (34) and (16), we shall refer to as the "equations 
of the gun." Equations (33) and (34) give 
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which, on treating /3 as a constant, yields 



m ''^ 



^-^Y^+cz{-E{A,V)-^kE}~\-A\ mx=0, (37) 

J-~P Jo 

using the conditions that when ^ = 0, V = 0, and E (A, P) = JcE. 

The equation (37) is the energy integral^ and may be used instead of (33) 
and (34) in the equations of the gun. 

We now again fall back on the provisional assumptions, and use the 
equation (30) and the law (24) for the rate of burning function. Equa- 
tions (34) and (16) may therefore be written, on treating k as constant, 

&n ~ = (C + Ax-^ ^4-7AP, (38) 

cix \ 0/ cm 

Clt Vq 

the initial conditions being that motion begins at a^ = when P = Eo, the 
value of E at a? — 0. Up to this stage we have 

m^ = (C/c-l/S)P. (40) 

When maximum pressure occurs during burning, we have from (38) 

^cTI dzjclx = ^AP, 

and therefore, from (32) and (39) 

-^ = Ml fo (,,, + , ,,/F), (41) 

on using accents to denote evaluations under the circumstances of maximum 

pressure. 

Dimensional Relations. 
9. Take 

^ = 1+^. (42) I' = c^s-^> (^5) 

Y = (l-.^)Al!L.v, (44) ,^ = 2z:^.^, (47) 



and for the time instead of t take 

% /C 1' 



t. (48) 



Then the equations of the gun may be written 

dd/cU = M.v, (49) dz/de = ddplcW + 'yp, (51) 

dv/dt = p — r, (50) dzjdt = {p-{--!^)^{z), (52) 
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the initial conditions being that 



^ = 1, 'y = 0, p •= z — 1% (53) 

while the energy integral becomes 

-~ {z-^^e) = J M^H ^rcie, (54) 

and the condition for maximum pressure when it occurs during burning 



The solution of the equations of the gun is then seen to depend on the 
constants M and tjy, and on the functions ^ (z) and r. These may remain the 
same over a considerable variety of circumstances from the point of view of 
applications, so that a gun calculation for one set of circumstances may be 
made to yield the solution for certain other sets of circumstances. This 
principle provides the proof of the general theorem immediately following. 

If A, C—- c/S, Xy m, e, Vq/sq, kU, ui, Uq, l — y8, 

be multiplied by a parameter raised respectively to the powers 

(A), (C), (x), (m), (c), (vo), (U), (m), («o), (/8), 

then V will be multiplied by the parameter raised to the power 

(A)-f(%) + (;8)~(m)-(^i), 

P will be multiplied by the parameter raised to the power 

(c)+(n)-(C), 

time will be multiplied by the parameter raised to the power 

W+(C)-(wi)-(n)-(c), 

providing (A)-{-(x) = (C), 

(c) + (m) + 2 (^1) + (H) = 2 (A) + 2 (vo) + (/3) 

(n) + (c) + (%) = (C) + (^o), 

and R be changed like P. 

Certain special cases of this theorem have long been used ; thus take 

(A) = 2, (C) = 3, (x) = l, (m) = 3, (c) = 3, (^0) = 1, 

(H) = (u,) = (^^o) = (/3) = 0, 

and we have the theorem that similar guns similarly charged with the same 
propellant produce equal pressures and equal velocities in linearly similar 
times. 
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General Differential Equation. 
10. The equations (49), (50), and (52) give 

a \ a^ / 1- rk \ 

and then on taking r constant during hand engraving, we obtain 

az X at/ . / . \ , / f r\\ 

On solving this equation for z we have, say, 

where the function Z depends on the function (p(z). Substituting in (54) we 
obtain, on rearrangement, 

^^ + 4(7-l)(f J -MZ |,.o, l^ + (t. + ro)^} +yMne = (y-l)Mn, 

(61) 
which is a differential equation for in terms of t. 

For pieces of propellant with a constant burning surf ace, so that <j^(2;) = 1, 
we have 

and the differential equation (61) becomes 

^^ + i(7-l)(^T-^ + 7Mro(^-l) = M(tiF + ro)^. (63) 

After band engraving, we take r = and obtain 

'^ (^0 1 cf^ 



i 



<l>{z) M dt 



+• w^ — (t;i + mti\ (64) 



where the suffix 1 refers to the end of band engraving. This gives 

z = Z ■^Zu^—+xnt — (vi + mti)i-, (65) 

and on substituting in (54) we obtain, on rearrangement, 

+ (7 - 1) Mro (6'i - 1) = 0. (66) 
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Again, for pieces of propellant with a constant Inrning surface, this 
differential equation becomes 

(67) 

When CT = so that the rate of burning of the propellant is proportional 
to the pressure, the method of the present section is not the most direct, and 
when CT 9^ so that the rate of burning is not proportional to the pressure, 
the differential equations we have obtained cannot be integrated in terms of 
elementary functions. Of course, in any particular case we may obtain a 
numerical solution to any required degree of accuracy, by the methods of 
step-by-step integration, and in this way tables could be constructed giving 
the solution in terms of M, xs and ro for each form of ^ {z). But the trans- 
formations of the present section make a contribution to the results of the 
next section. 

Concentrated Band, Characteristic Parameters. 

11. By concentrated driving band we mean the conditions presented by 

r = 0, for ^ > 1 ; ro 9^ 0, 

which is Charbonnier's method of treating band engraving. 

The solution of the equations of the gun is then determined completely by 
the function <f> {z) and the values of 

M, vs, ro, 6^3, (68) 

the suffix 3 referring to the muzzle. In particular, the value of v^ and 
whether the charge is all burnt or not, depend only on these. If the charge 
is all burnt in the gun, then whether the maximum pressure occurs while it 
is burning or not, depends only on ^{z) and the values of 

M, m, ro, (69) 

as does also the value of p\ 

For any one form of the function <p (z) the functional relations here 
indicated are the same for all guns and propellants to which the equations 
(33), (38), (39) apply. 

When (^ {z) = 1, and maximum pressure occurs during burning, we see from 
(67) that Mp' depends only on Mot and Mro, since for a concentrated band 

01 = 1, vi = 0, zi=z ro. 

Eeverting now to the notation of § 8, we see that the solution of the 
equations of the gun, and in particular the value of 

(l--y8)A^0 ^ ^ 
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depend only on 

(l-^)A^^o^ ^ {G-^clh)no ^ ^:=^Eo, C3-^ /g /7i^ 

C3 denoting the total capacity of the gun. 

For the same propellant whose rate of burning is not proportional to the 
pressure, these may be replaced respectively by 



mso 



(1 — y8)A'yo 



' . V3, (70)' 



and ^^^=^^. ^, Eo, ?F^. (71)' 

When the charge is all burnt in the gun, then whether maximum pressure 
occurs during burning or not, and also the value of 

^ F. (72) 

depend only on the values of 

(l~;g)A^V (C-£/S)2^o^ ^Zl^Eo. (73) 

Again, for the same propellant whose rate of burning is not proportional to 
the pressure, these may be replaced by 

P' (72)' 

and a^mp^^ C-^, B„, (73)' 

kcmS(f kc 

respectively. 

If, further, the pieces of propellant have a constant burning surface, and 
maximum pressure occurs during burning, then 

(l-/Q)(C-o/8) / A% y p, ,^ 

m xkcliuisj 

depends only on 

7?i, \kcJliiisJ Ui m \kcJluisJ 

One of the relations of this section has been used by F. B. Pidduck. Taking 
the same propellant with the rate of burning proportional to the pressure, 
.so that 2/.0 = 0, we have that 

(Cz:f/i) p', (72)'' 

roG 

depends only on 

(1-/3) A^^o^ (C-~c/8) ,, .^..„ 

kcms(f kc 



The Principles of Internal Ballistics, 



303 



Pidduck plotted observed values of the first of these three quantities 
against the second, and sought the effect of the third. 

The method appears capable of much further utility. For example, on 
taking the same propellant and 

^{Z) = 1, ItQ = 0, 

we have, when maximum pressure occurs during burning 



depending only on 



m \kcsj 

(l-/ 3)(C-c/S) /A^oY ^^ 



(74)' 



(75)' 



m \kcso/ 

But the greatest value of the present relations is when the rate of burning 

is not proportional to the pressure, and the equations of the gun cannot be 

integrated by elementary functions. For example, when with the same 

propellant 

(j) (z) = 1, uo=^ 0, 

and when the maximum pressure occurs during burning, P' depends only on 

(l^/3)(C^c/5) ^AV 
m 



\2 
\kG So/ 



(75/' 



After Engraving wit/i Bate of Burning Proportional to Pressure. 
12. We shall now take r = 0, ra- = 0, and start from the conditions 

= 0u Z=: Zi, V = Vi, i? =i?i. 

From (50) and (52) we have 

dv/dz = l/^(z), 

V = Vi-{- \ 



so that 

and from (49) and (52) we have 



Mi; dz 



<i> {z) dd' 
On substituting in (54) from (77) and (78) we obtain 

' dz 



(76) 

(77) 

(78) 



4(^-i)|S = 



Vi-\' 



j 



^(2) 



2_ 

M 



{;^-|M-{''.+ 



dz 



2l 



4>{z). 



H^r 



giving 







-^1 + 



dz 



r, 4> (^) 



z 



M 1 7-1 



rdO y -- \ Vi-\- 



r 



dz 



2 



dz 
J{z) 



. (79) 



<^(^) 



which is the general 0, z relation after engraving and during burning. 
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For a concentrated band (77) and (79) reduce to 

*^ dz 



V 



<"' clz 



(77)' 



and 



log^^(Y~i) 



r, 



■o^(^) 



2z 



dz 



dz 

2 • ^(zi 



(79)' 



,,(7-l)M Ur„</.(«), 

The equations (77)' and (79)' form the basis of Charbonnier's solution, but 
in evaluating (79)' he approximates except for <^ {z) = 1. 
On evaluating (79) for 

<i>{z)=l, ^(z)={l-z)i, 



we obtain 










(80) 



with 



a = (b-j-c^)^, 



where for ^ (;<;) = 1, or a constant burning surface, 



^=^> ?i 



^1, 



^ 



2_ 
M 



■^1 
1 



e=l/{j^l)M^(zi^v,) ^' 
X = 1/(7-1) 

and for <^ {z) = (\,—%f, or cylindrical pieces of propellant, 

K^{\-z)\ ^i = (i_«i)i 

1 r 1 / 1 



(81) 



■^ 



Z^ = 



C 



\ 



/^ 



l-fl/2(7-l)Ml2M\7~-l 

g'l + l^i 
l + l/2(7-l)M 

1 
7-1 + 1/2M 

2M ?i + -|^'i 



r^^ -(fi + J^i)3 



^. 



(82) 



^ 



c^ {l + 2(7-l)Mp 

After all the charge is burnt, our equations reduce to those in general use. 
The result for a concentrated band may be obtained very easily by putting 

= 1, % = 0, ^1 = ro {^ 0), 

in (81) and (82), but the derivation of those for no hand requires a limiting 
process, so that they are most easily obtained directly. 
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For <^(^) = 1 the results are given by Oharbonnier,* while for ^(z) = (1— «)i 

we have 

v=2{l^(l^z)i}, (83) 

pd = [^-^2(r^-^l)M {l-(l-0)i}]. (85) 



The Harmonic Development of the Tide -generating Potential 
By A. T. DooDSOK, D.Sc, Tidal Institute, the University of Liverpool 

(Communicated by Prof. H. Lamb, F.E.S. Eeceived June 12, 1921.) 

§ 1. Introduetion. 

The harmonic development of the tide-generating potential is the basis of 
most work on tidal observations, and since 1883 the development given by 
Sir G. H. Darwin has been universally used and has been of remarkable value. 
But discrepancies between prediction and observation are serious and have 
been attributed to faulty " harmonic constants " ; it has been assumed that if 
these were improved better predictions would be obtained, and it has also 
been tacitly assumed that it is only necessary to consider the harmonic 
constituents as given by Darwin. Eecent work, however, especially at the 
Tidal Institute, has shown that when all the '' Darwinian constituents '' are 
removed from the tidal height there is a residue composed of constituents 
which are not included in his schedules. These are such that any slight 
improvements possible in the " constants " usually obtained are comparatively 
negligible. 

The obvious course, therefore, was to make a more thorough development 
of the potential, and in view of the unknown nature of the residues great 
accuracy was obviously desirable, especially as the possibility of resonance has 
always to be considered. The development given in this paper, even if it 
proves to be needlessly thorough for practical tidal work, will cover the needs 
of research work, since it includes all terms whose coefficients (relatively to 
the greatest coefficient) are greater than 0*00010. 

Darwin used the old lunar theory and referred everything to the orbit 
rather than to the ecliptic ; his results are all given in the algebraic form, 
arithmetic being used only to decide what terms to omit. His development 
is a quasi-harmonic development because he retains factors in the coefficients 

^ Loe, city Chap. VII. 
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